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Abstract 

We study the properties of elongated Bose-Einstein condensates. First, we show that the 
, dimensions of the condensate after expansion differs from the 3D Thomas-Fermi regime. 

We also study the coherence length of such elongated condensates. 

1 INTRODUCTION 

, The experimental study of one dimensional (ID) degenerate Bose gases, where the radial motion 
■ is "frozen", is currently an important direction of research in ultra-cold atom physics |1I2||^ . In 
uniform ID systems, fluctuations of the condensate phase are pronounced, because of a large 
population of low-lying states; as a result, no long range order exists at any temperature. In a 
. trapped ID gas, the finite size of the sample naturally introduces a low-energy cutoff, and at 
' sufficiently low temperature a phase coherent sample can exist (see W and references therein 
. for the properties of ID gases). Essentially, the same analysis holds for three dimensional (3D) 
' trapped gases in very elongated traps |S] : here the interaction energy per particle is larger than 
the radial frequency, and the condensate wavefunction is built from several radial modes. Still, 
low-energy excitations with a frequency smaller than the radial frequency are effectively one- 
dimensional, and induce large fluctuations of the condensate phase |S| . Such condensates with 
a fluctuating phase, or quasicondensates, have been observed experimentally in |Sl[7ll5l [Tni m| . 

We present here an overview of our experimental study of elongated condensates in a very 
anisotropic trap. Our set-up, presented in section |2] allows us to characterize precisely the 
properties of degenerate gases between the 3D and ID regimes. In this regime, we first show 
(section|21l that the (quasi)condensates exhibit the same smooth profile as a true (fully coherent) 
condensate. This profile, however, is affected by both the radial quantum pressure (zero point 
oscillation in the tight trapping potential) and the presence of an interacting thermal cloud 
surrounding the condensate and interacting with it through two-body repulsive interactions. 
In section ^ we discuss the coherence length of such elongated condensates. Using Bragg 
spectroscopy |13lll0lll2llllj . we show in section ETI that the coherence length is indeed strongly 
affected by large fluctuations of the condensate phase and that the measured momentum 
width depends strongly on temperature; the corresponding coherence length L^, in the ratio 
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5 < ^ < 30, is in good agreement with theory [Sj^]. For larger coherence length, Bragg 
spectroscopy becomes increasingly difficult to apply, and a second method, using a 7r/2 — 7r/2 
interferometer similar to [S] is presented in section [4. 21 



2 CREATION OF ELONGATED BEC IN A IRON-CORE ELECTROMAGNET 

The experimental setup has been described elsewhere Briefly, a Zeeman-slowed atomic 

beam of ^^Rb is trapped in a MOT, and after optical pumping into the F — 1 state is transferred 
to a magnetic loffe-Pritchard trap created by an iron-core electromagnet. A new design allows 
us to lower the bias field to a few Gauss and thus to obtain very tight radial confinement jJS] . 
Final radial and axial trap frequencies {ujj_ and lOz) are respectively 760 Hz and 5 Hz (i.e. 
aspect ratio = 152). In this way, we obtain needle-shaped condensates containing around 
5 X 10^ atoms with a typical half-length L ~ 130 fim and radius R ~ 0.8 ^m. Since 
the chemical potentiel ^ 4 (see definition in section |2Jl, the clouds are between the 3D 
and ID Thomas- Fermi (TF) regimes Jji'i2ii^- The measurement of the number of atoms is 
calibrated from a measurement of Tc ■ The temperature T is extracted from a two-component 
fit to the absorption images or is extrapolated from the final frequency of the rf ramp when the 
thermal fraction is indiscernible. After the last evaporation ramp, we observe the formation of 
a condensate that exhibits strong shape oscillations (see jT2] and references therein) . This is in 
line with a recent observation by Shvarchuck et al. ^ , who propose the picture of a BEC locally 
in equilibrium at each point on the axis, but globally out of equilibrium. Indeed, the radial 
degrees of freedom are rapidly thermalized, but axial motion damps on a much longer time 
scale. We therefore hold the condensate for a variable time (typically 7 seconds) in presence 
of an rf shield, so that these collective oscillations have no influence on the Bragg spectra. 



3 EXPANSION OF ELONGATED CONDENSATES 



In anisotropic traps, the aspect ratio of the ballistically expanding cloud can be used to find 
a signature of Bose-Einstein Condensation. For a Bose gas above critical temperature, the 
expanding cloud is expected to become isotropic at long expansion times (typically for times 
larger than where ojz is the smaller trap oscillation frequency). On the contrary, a Bose- 

Einstein condensate of initial aspect ratio e = liJz/uj± will only become isotropic at a fixed time 
~ wj^ and then will invert its anisotropy for longer expansion time. 

The transition from a 3D to a lower dimension (ID) condensate is also observable through 
the expansion of the degenerate gas T. For a Bose-Einstein condensate in the 3D Thomas- 
Fermi regime (/i ^ ?iwj_), the half-length LTF(long axis) and the radius i?TF (confined di- 
mension) of the condensate vary similarly with the number of atoms iVo. More precisely, one 
finds the 3D Thomas- Fermi length L^p = 2fj,/{Muj'^) and radius Rtf = c^tf, with n given 
by 2^i = to(15iVoa/cr)2/5 where Q = (luIlUz)'^/^ , a = ^/h/{MiD) and a = 5.32 nm jH]. As 
a result, the initial degree of anisotropy e = R/L = lOz/loi_ ^ 1 is independent of the num- 
ber of atoms. After a period of free expansion of duration t, the dimension of the expanding 
condensate are found using the scaling solutions jl9j: 

i?(t) =i?TFVTT^. (3.1) 



i(f) = LTF(l + e^ 



r arctan r — In yl 



TF- 



(3.2) 



with the dimensionlcss time variable t = uj±t. Equations H3. 113.2(1 show that the aspect ratio 
is independent of the atom number at any time, and approaches Lot for r ^ 1 
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Fig. 1. Half-length L (upper-left), radius R (lower-left) and aspect ration (right) of the quasi-lD 
condensate versus the condensed atom number after a 24.3 ms time of flight. Experimental values 
(diamonds) are compared to the standard T = Thomas- Fermi (plain), and the self consistent solution 
of equations 13.41 and 13.31 (dashed). The effect of quantum pressure alone is also shown (dot-dashed). 
Calculations use the measured temperature and number of condensed atoms, as determined by a two- 
component fit to averaged absorption images. Combined statistical and systematic errors are indicated 
with the confidence band (grayed area) around the experimental points. 

As one pushes towards the ID regime (/z ~ huj±), the radial expansion is increasingly 
affected by radial zero-point motion in the trap, neglected in the 3D TF approximation. Ulti- 
mately, the radial size in the trap should approach a± — y^h/Mu!±, and the radial expansion 
velocity uj±a± = ^Jhjj^^jM. Since axially the TF approximation (/x ^ fiwz) is still valid, the 
length still depends on A^o, and so does the aspect ratio. Accordingly, a dependence of the 
aspect ratio on A'o at near-zero temperature can be seen as a signature of the breakdown of 
the TF approximation in the radial direction 

To show how quantum pressure manifests itself in our experiment, we first turn to the 
measurement of the axial half-length L after release, which faithfully reflects the half-length in 
the trap since axial expansion is negligible (eq. I3.2|l for our experimental parameter. Figure ^ 
shows that the measured values of L for various number of condensed atoms (filled diamonds) 
are smaller than the standard Thomas- Fermi prediction (solid line) . This smaller value results 
from radial quantum pressure (ID), and to a slightly larger extent from the compression of 
the quasi-condensate by the 3D thermal component (excitations with energy much larger than 

We take into account radial quantum pressure in this crossover regime by finding the optimal 
ground state energy 

I [E^V7l+^(g.W/^n(v/T^-0)^/^]. (3.3) 

i—\ i—1 

where are effective parameters taking the zero-point energy contribution in the direction 
i = {x,y,z} into account. Solving (eq. 13. 3|) allows then to calculate a modified aspect ratio 
in the trap e — jz/jx,y- To include the effect of the thermal cloud, we use the Hartree-Fock 
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approach [2011 17| and calculate the modified length : 



1(0) 



(3.4) 



with the coupling constant g = 47r?i^a/M, the thermal de Broglie wavelength At = [27rS^/(AffcBr)]^/^, 
and g3/2(;c) = x^/n'^^^- The black dashed fine in Fig.^shows the solution of Eq. (|3.4(l 

assuming a paraboHc profile, such that no(0) = 15NoL~^e^ / {Sir). The modified aspect ratio in 
the trap e is calculated according to l3.3l and by a self-consistent iteration of equations H3.4(l and 
we obtain L and "ix,y,z- We find that the calculated lengths are close to our measurements 
to within our estimated calibration uncertainty (shown by the grayed area explained in the 
caption). 

We then turn to the measurement of the cloud radius after 24.3 ms of time of flight. Figure^ 
shows that the measured values of R at various condensed atom number (filled diamonds) are 
larger by about 20% than the standard Thomas-Fermi prediction (filled line). This increased 
radius results again from quantum pressure, which leads to an excess released energy in the 
transverse direction. Using again [21], the radius after time of flight is found to be: 

with 7i calculated from (eq. I3.3|l . Calculating the effect of the compression of the thermal cloud 
and the subsequent evolution in time of flight is more complex, since it involves energy exchange 
between both components. By making the extreme assumption that all the interaction energy 
is converted to kinetic energy into the condensate after expansion, we can use again eauation l3.5l 
with correction from 13.41 and 13. 31 The results presented in figure ^ show that, both approaches 
agree with the experiment to within our calibration errors. 

Finally, we observe a clear dependence of the aspect ratio on the condensed atom number, 
which can be correctly reproduced by our theoretical approach. This proves that in our trap, 
with this atom number ('--^ 5 x 10''), the condensate is at the onset of ID. 



4 MEASUREMENT OF THE ID PHASE FLUCTUATIONS 



Phase coherence of the (quasi)condensate can be characterized by the 1/e decay length of the 
spatial correlation function, C^^^\s) = 1/iV / d3R(*t(R + se,/2)\E'(R - se,/2)) E7|, with ^ 
the atomic field operator and the unit vector in the axial direction. Consider a sample 
of No condensed atoms, trapped in a cylindrically symmetric, harmonic trap with an aspect 
ratio A = uj^/ui±_ <C 1, in the 3D Thomas-Fermi regime. Petrov et al. investigated 
the long wavelength fluctuations of the phase of 4", with the following conclusions. Below a 
characteristic temperature = 157Vo(Sa'z)^/32/j,, the phase fluctuations are small, and the 
coherence length is essentially the condensate size L, as in "ordinary" 3D condensates 13',!22| 
(to put it differently, C^^^ is limited by the density envelope of 5*). On the other hand, for 
T ^ the coherence length is limited to a value much smaller than L, the value at the center 
of the trap being = LT^/T. 

The theoretical results obtained in 5,,l2 do not apply directly to our experiment. In those 
works, a 3D Thomas-Fermi density profile was assumed, whereas in our case modification of 
the density profile by the thermal cloud and radial quantum pressure must be accounted for. 
As discussed in the previous section, we indeed find that a parabolic profile is still a good fit 
function, but that the usual T = relations between /i, L, R, and the number of condensed 
atoms iVo are no longer valid. We therefore extend the calculation of the axial correlation 
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function in jl2| to an arbitrary density profile, in the local density approximation. In the mean 
field regime 111, from the result of for a ID uniform Bose gas at finite temperature, we 
obtain: 

C(.)=/..n.(.)exp(^£iM-), (4.1) 

where ni{z) — J (Pr±no{r±, z) is the axial ID density of the quasi-condensate, while no{r) 
is its 3D density profile. The coherence length near the center of the trap is given by — 
h^ni{0)/{MkBT) Following Petrov et al. |4I5 | . we define the temperature which delineates 
the border between coherent and phase-fluctuating condensates as = L^T/L. Since ni(0), 
L, and T are extracted directly from the images, the definitions of and relate the 
coherence properties to experimentally measured quantities. The axial momentum distribution 
follows from a Fourier transform of C(s) and is well approximated by a Lorentzian of width 
Ap0 = ah/L^ (HWHM), with a = 0.67 for a 3D Thomas- Fermi condensate [Hll^ . 



4.1 Momentum spectroscopy with Bragg difFraction 

A convenient technique to measure the coherence length is to analyse the momentum distri- 
bution along the long direction of the condensate. Our measurement is based on four-photon 
velocity-selective Bragg diffraction. Atoms are diffracted out of the condensate by interaction 
with a moving standing wave, formed by two counter-propagating laser beams with a relative 
detuning S T3 26 . Due to the Doppler effect, the momentum component resonantly diffracted 
out of the condensate is Pz = M{5 — 8ajR)/(2A:L) with uj-r. = f^k^/ i'^M), M the atomic mass, 
and /cl = 27r/A (A = 780.2 nm). The predicted spectral width is therefore aAj/^, where 

2/ifcL 



To build the momentum spectrum of the quasi-condensate, we measure the fraction of 
diffracted atoms versus the detuning 5 between the counter-propagating laser beams. The 
differential frequency 5 must be stable to better than the desired spectral resolution, about 
200 Hz for our typical L^p — 10 /im. The optical setup is as follows. A single laser beam 
is spatially filtered by a fiber optic, separated into two arms with orthogonal polarizations, 
frequency shifted by two independent 80 MHz acousto-optic modulators, and recombined. The 
modulators are driven by two synthesizers stable to better than 1 Hz over the typical acquisition 
time of a spectrum. The overlapping, recombined beams are then sent through the vacuum 
cell, parallel (to within 1 mrad) to the long axis of the trap, and retro-reflected to obtain two 
standing waves with orthogonal polarizations, moving in opposite directions. The lasers are 
tuned 6.6 GHz below resonance to avoid Rayleigh scattering. The laser intensities (about 2 
mW/cm^) are adjusted to keep the diffraction efficiency below 20%. 

Bragg spectra have been taken at various temperatures between 90(10) nK and 350(20) nK, 
while Tc varied from 280(15) nK to 380(20) nK. The temperature was fixed to within 20 nK by 
controlling the final trap depth to a precision of 2 kHz. The line shape of the spectra is clearly 
Lorentzian, not Gaussian (see Fig. 12). This is a significant result, because a Lorentzian-like 
profile is expected for a momentum distribution dominated by phase fluctuations (see ^2]), in 
contrast to the gaussian-like proflle expected for a pure condensate |27I13| . From the Lorentzian 
flt, we extract the measured half-width Ai^m for each temperature. 

Figure 121 shows the measured spectral width A^m versus Ai^^. Dispersion of the measured 
ni(0) and T results in the horizontal error bars on Aj/^, while vertical error bars indicate the 
standard deviation of the fit width. We assume that all experimental broadenings result in a 
Gaussian apparatus function of half- width wq , to be convolved by the Lorentzian momentum 
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Fig. 2. (a) Absorption image of degenerate cloud (center) and diffracted atoms (left and right), 
averaged over several shots, after free-flight expansion, (b) Diffraction efficiency (rescaled) versus 
relative detuning of the Bragg lasers at T = 261(13) nK, corresponding to T /T^ — 20(2). A typical 
statistical error bar is shown. This spectrum is the superposition of 12 "elementary spectra", as 
described in the text. The true average center is 30.18(2) kHz, close to 30.17 kHz, the four-photon 
recoil frequency. The solid line is a Lorentzian fit, giving a half-width at half-maximum (HWHM) of 
316(10) Hz. (c) and (d) show respectively the residual of a Lorentzian and of a Gaussian fit to the 
above spectrum. Residuals are folded around the 5 = axis, and smoothed with a six-point-wide 
sliding average. 
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Fig. 3. Left : Half-widths at half-maximum Ai/m of the experimental Bragg spectra versus the pa- 
rameter Af^ oc /i/I/0 (see text). Vertical error bars are the standard deviations of the fit width; the 
horizontal error bars are the one-sigma statistical dispersions of Ai/,/,. The solid line is a fit assuming 
a Voigt profile for the spectra. Right : coherence length deduced from the width, after correction for 
the gaussian broadening, 



profile with a half- width aAv^j,. The convolution, a Voigt profile, has a half-width oiAv^j^ -|- 
a/uiq -t- (aAv^y- 1 \. Note that fitting a Voigt profile instead of a Lorentzian to a spectrum 
gives the same total HWHM to less than 5 %, but the Lorentzian shape is too predominant to 
extract reliably the Gaussian and the Lorentzian contributions to the profile. Using a and 
as free parameters to fit the data of Fig. 13 we find wq — 176(6) Hz, and a = 0.64(5)(5). The 
first uncertainty quoted for a is the standard deviation of the fit value. The second results 
from calibration uncertainties on the magnification of the imaging system and on the total 
atom number, which do not affect wq- The agreement of the measured value of a with the 
theoretical value 0.67, to within the 15% experimental uncertainty, confirms quantitatively 
the temperature dependence of the momentum width predicted in Ref. j^]. The coherence 
length hj-p^ deduced from this measurement varies between 5.9(8) and 39(4) /zm, in the range 
6 < TjT^ < 28. 
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Fig. 4. Left : The atom interferometer used for our coherence length measurement uses 2 7r/2 pulse to 
create and overlap 2 displaced copies of the same condensate. Right : The fringe spacing and contrast 
changes as the time between the 2 pulses (i.e. the separation between the 2 copies) is changed. A 
measurement of the fringe contrast directly measures the spatial first-order correlation function of the 
matter-wave field. 

4.2 Interferometry with elongated BECs 

The Bragg spectroscopy method has proven to be appropriate for measuring short coherence 
lengths. However, when applied to nearly coherent condensates {Lc L), the extremely nar- 
row Bragg resonance becomes increasingly difficult to resolve (for L$ ~ 40/im, the momentum 
width is already of the order of our resolution). However, since our goal is to investigate the 
cross-over from a phase coherent to a phase-fluctuating condensate, we need to measure co- 
herence length larger than the 39(4) /im attainable with Bragg spectroscopy. For this, we have 
implemented a matter- wave interferometer ^ . As shown on figure ^ this interferometer uses 
a first ^ bragg pulse to split the condensate in 2 equal copies that fly apart from each other 
at a velocity 2vji = 2-|y ~ 1.2 cm/s. After a time of flight 2ms < T < 10ms, the 2 copies, that 
separate by 2vii x T, are recombined using a second ^ pulse, that creates 2 interferometer out- 
put (see flg. For finite T, complementary fringes appear in the two output ports, because 
of the smoothly varying expansion phase that results from the conversion of a small fraction 
of the initial mean-field energy e^/z) into axial kinetic energy. The resulting fringe pattern 
is Fourier analyzed to find the contrast and the fringe spacing in reciprocal space. In this way, 
we filter out shot-to-shot, global phase fiuctuations that would result in a translation of the 
fringe pattern as a whole. This method allows for averaging many processed images without 
artificially reducing the contrast. After this processing, this method allows a direct measure- 




separation ( m) Separation ( m) 



Fig. 5. Example of 2 direct measurements of the correlation function (fringe contrast) as a function 
of the separation s between the two copies of the condensate ; left ; Tg>/T — 0.45 and 2 x 10^ atoms 
(Ltp ~ lOO^m), right : T^/T = 0.75 and 1.3 x 10^ atoms (Ltp ~ 90/xm) 
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ment of the correlation function C'^^-' as a function of separation s (see figure^, irrespective of 
the amount of phase fluctuations. 

5 CONCLUSION 

We have demonstrated two important features of elongated Bose-Einstein condensates: (i) the 
dimension of the condensate after expansion; (ii) the momentum distribution of quasiconden- 
sates; Our method, along with atom interferometry can be applied to investigate how long 
range order develops during the condensate growth. 
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